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A method for k-space data acquisition and MRI device 




FIELD OF THE INVENTION 

The present invention is related to the field of magnetic resonance imaging 
(MRI), and more particularly, to k-space data acquisition, 

• • * 

5 BACKGROUND AND PRIOR ART 



between spatial encoding (phase encoding and firequency encoding in the time domain) and 
the corresponding image obtained by applying the Fourier transform. For data acquisition in 
the k-space a sampling trajcctoiy of a frequency-encoded signal is typically used. The basic 
10 concepts of the k-space notation are explained in more detail in "Principles of Magnetic 
Resonance Imaging, a Signal Processing Perspective", Zhi-Pei Liang, Paul C. Lauterbur, 
IEEE Press Series in Biomedical Engineering, 2000, in particular chapter 5.2.3, pp. 157. 



J. Phys. C: Solid State Phys,, vol. 10, pp. L55-L58, 1977 an MRI method is known which is 
15 commonly referred to as echo-planar imaging (EPI). The term EPI is broadly used to refer to 
the class of high-speed imaging methods that collect a "complete" set of two dimensional 
encodings during the free induction decay period following a single excitation pulse. Hence, 
EPI has become a synonym for single-shot imaging, although multi shot EPI methods with 
interlaced k-space coverage are also in common use. 
20 A key concept of EPI is the use of time- varying gradients to traverse k-space. 

For k-space data acquisition a variety of trajectories are known fi'om the prior art such as zig- 
zag trajectory, rectilinear trajectory and spiral trajectory. For a discussion of the various prior 
art trajectories reference is made to the above referenced book of Liang and Lauterbur, 
chapter 93, pp. 303. 

25 Figure 1 and Figiae 2 show rectilinear trajectories for k-space EPI data 

acquisition. 

In a first data acquisition step a first partial data acquisition of the k-space of 
the target region is obtained by following the trajectory as depicted in Figure 1 . The 
trajectory starts at the central point 100 of the k-space. From ttiere it goes into the lower left 



The k-space notation is widely used in the art of MRI to establish a connection 



From P. Mansfield, "Multi-planar image formation using NMR spin echoes,' 
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comer 102 of the k-space. Starting fiom the lower left comer the k-space is partially scanned 
by means of a rectilinear trajectory. As it is known from WO 02/35 250 Al half of the k- 
space is covered plus an additional seven lines. 

After the first partial k-space acquisition in accordance with Figure 1 a brief z- 
5 shimming gradient pulse z - zo is applied before the second partial k-space data acquisition is 
performed in accordance with the trajectory of Figure 2. 

TGe^^eclofy "of FiguiF2lTalsoTectilinearm^ 
where the trajectory of the partial k-space acquisition of Figure 1 ends. 

• A k-space acquisition scheme of the type shown in Figures 1 and 2 is known 
10 from WO 02 / 35 250 Al and from "Single-Shot EPI With Signal Recovery From the 

Susceptibility-Induced Losses", Allen W. Song, Magnetic Resonance in Medicine 46:407- 
41 1 (2001) for application in functional magnetic resonance imaging (fMRI). Based on each 
of the partial k-space data acquistions an image is obtained and the two images are 
overlapped in order to produce a resulting image. 
15 It is an object of the present invention to provide for an improved method for 

k-space data acquisition in order to increase the spatial sensitivity of MRI. 




SUMMARY OF THE INVENTION 

The present invention provides a method, a MRI device and a computer 
20 program product featuring improved k-space data acquisition for increased spatial sensitivity 
of the acquired images. 

The present invention is based on the discovery that a forward and subsequent 
reverse k-space acquisition (cf Figure 1 and Figure 2) leads to brain activations as detected 
with fMRI in very different places in the target region. In particular this decreases the 
localisation and sensitivity of the resulting images in the blood oxygenation level-dependent 
(BOLD) statistical analysis. 

The present invention provides a method for k-space data acquisition which 
overcomes this disadvantage of the prior art. In essence, the data acquisition in at least two k- 
space is performed in parallel in an interleaved manner. For example after sampling into a 
first direction in the first k-space a compensation pulse is applied. A sampling is performed in 
the second k-space in the opposite direction. In the next step another compensation pulse is 
applied before the first k-space is again sampled in the first direction. The at least two k- 
spaces cover the same physical region, such as a slice of a patient's body. Individual images, 
each having its own information content, are associated with their respective k-spaces. 
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Within the framework of the present application the term compensation pulse 
indicates a pulse which affects artefacts in the reconstructed image. In particular such 
compensation pulses affect susceptibility artefacts. Very good results as to compensation of 
susceptibility artefacts artefacts. Very good results as to compensation of susceptibility 
5 artefacts are obtained when magnetic gradient pulses, such as 2-sIdnned gradient pulses are 
employed as compensation pulses. 

In accordance with a preferred embodiment of the invention an image is 
formed based on the data samples which have been acquired for each of the at least two k- 
spaces. The images are combined to fonn a resulting image whirfi features improved 
10 localisation and sensitivity. 

The invention is particularly advantageoxis in that the echo times (TEs) for the 
two or more k-spaces do only have a small difference of for example less than 1 msec. This 
compares to the prior art, especially the Song reference, where even in the case of partial k- 
space coverage of the two echoes, the TEs differ by some eight to ten milliseconds. This 
15 results in an amplitude difference of the echo signals which affects the quality of the 

composite image, for both Sum-of SQuares (SSQ) or Maximum Intensity Projection (MSP) 
combination. Typically, T2* (transverse relaxation time) being approximately 50 and 30 ms 
for main magnetic field strengths of 1 ,5T and 3T, respectively, the amplitude difference 
between the two images at an echo time difference of 10 ms will amount up to 30-40%, while 
20 BOLD related signal variations are typically on the order of 5-10%. Even amplitude 

correction will hamper adequate composite image formation in the prior art due to the image 
contrast mechanism of the BOLD effect. 

As a consequence the small difference in the TEs which is accomplished in 
accordance with the invention has the effect to greatly enhance the quality of the composite 
25 image. This makes Uie invention particularly suitable for functional magnetic resonance 
imaging (fMRI). 

It is a further advantage of the present invention that typical BPI corrections to 
line up echoes from positive and negative read out gradients do not need to be performed as 
all echoes contributing to one image have the same polarity. This improves the robustness of 
30 the IQ per "echo" image. 

In accordance with a further preferred embodiment of the invention a larger 
number of k-spaces is utilized for the data acquisition. An even number of k-spaces is 
preferred as this has the advantage of not having to use the EPl phase correction when 
acquiring in both directions. Further it is preferred to apply a ky gradient pulse after each 
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horizontal scan or after every n-th scan when a number of n of k-spaces with n different 
compensation values is used. 

In accordance with a further pref^ed embodiment of the invention the TE 
difference for the two k-space data acquisitions is just the time of one horizontal readout 
5 along the kx axis. 

As any ky line for echo 1 is acquired immediately (one echo spacing, typically 

1):5— 0:8Tns)T)riortothe-samerline-forccho^rtheTesuk 

images'* is just a single echo spacing apart. Consequently, tiie contrast (amplitude) of the 
images is the same, and combined usage of SSQ or MIP works much better in comparison to 
10 the Song method. 

In accordance with a further preferred embodiment of the invention the echo 
images are acquired with a partial Fourier coverage methodology. In particular the prior art 
problems with differing flow sensitivities in backward and forward partial Fourier EPI 
readouts can be resolved by application of the disclosed method. 
15 In accordance with a further preferred embodiment of the invention the 

starting point for the sampling is the center of the k-spacc. From there the trajectory goes to 
the lower left comer of the first k-space. The following scans in the first k~space always go 
into the same direction. After each sampling along the kx axis of the first k-space a 
compensation pulse such as a z-shimming pulse is applied before a sampling into the 
20 opposite direction is performed for the second k-space. This way data acquisition for the two 
k-spaces is performed concurrently in an interleaved manner. 

In accordance with a further preferred embodiment of the invention the ky 
position of the sampling is incremented after each sampling along the horizontal axis in k.^ or 
in -kx direction. Alternatively the ky position is incremented only after sampling into the k^ or 
25 -kx direction. In any case a compensation pulse is applied after each horizontal sampling in 
the kx and -kx directions for the interleaved data acquisition. 

In accordance with a further preferred embodiment of the invention the k- 
space data acquisition is performed partially for both k-spaces. This means that the combined 
trajectories of the k-space data acquisition cover half the k-space plus an additional couple of 
30 lines. This way the read out time is reduced. 

A further substantial advantage of the proposed method is that it enables to 
solve the problem of combining backward and forward partial fourier. This allows a shorter 
echo time, and it makes the functional contrast equal, especially by removing the difference 
in flow sensitivity. 
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In accordance with a further preferred embodiment of the invention the 
method of the invention is implemented in a computer program product. The computer 
program product can be used for the control unit of an MRI device in order to perform a 
method of the invention. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the following preferred embodiments of the invention will be described in 
greater detail by making reference to the drawings in which; 

J Figure 1 is a prior art example for partiaLk-space data acquisition, ► 

10 Figure 2 is a prior art example for partial k-space data acquisition which is 

subsequent to the data acquisition of Figure 1 after application of a compensation pulse, 
Figure 3 shows a MR devices having a control unit being programmed in 
accordance with an embodiment of a method of the invention. 

Figure 4 is illustrative of an embodiment of a k-space data acquisition method 
15 in accordance with the invention. 

Figure 5 is illustrative of a gradient switchmg diagram for the k-space data 
acquisition of Figure 4, 

Figure 6 is illustrative of a flowchart for k-space data acquisition in 
accordance with the methods of Figures 4 or 5, 
20 Figure 7 is illustrative of an embodiment of a k-space data acquisition method 

where a ky gradient pulse is applied after every second horizontal scan. 

Figure 8 is illustrative of a gradient switching diagram in accordance with the 
embodiment of Figure 7 of the inventive method, 

Figure 9 is illustrative of an embodiment of a k-space data acquisition method 
25 with a nimiber of n = 4 k-spaces. 

Figure 10 is illustrative of a gradient switching diagram in accordance with the 
k-space data acquisition method of Figure 9. 

DETAILED DESCRIPTION 
30 Figure 3 shows a magnetic resonance device 1 which includes a first magnet 

system 2 for generating a steady magnetic field, and also several gradient coils 3 for 
generating additional magnetic fields having a gradient in the X, Y, Z directions. The Z 
direction of the co-ordinate system shown corresponds to the direction of the steady magnetic 
field in the magnet system 2 by convention. The measuring co-ordinate system x, y, z to be 
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used can be chosen independently of the X, Y, Z system shown in FIG. 3. The gradient coils 
are fed by a power supply unit 4. An RF transmitter coil 5 serves to generate RF magnetic 
fields and is connected to an RF transmitter and modulator 6. 

A receiver coil is used to receive the magnetic resonance signal generated by 
the RF field in the object 7 to be examined, for example a human or animal body. This coil 
may be the same coil as the RF transmitter coil 5. Furthermore, the magnet system 2 encloses 

-arrexaminatiDn-space-whrchrislarge^ntiugh-^^ 

examined. The RF coil 5 is arranged around or on the part of the body 7 to be examined in 
this examination space. The RF transmitter ijoil 5 is connected to a signal amplifier and 
demodulation unit 10 via a transmission/reception circuit 9. 

The control unit 1 1 controls the RF transmitter and modulator 6 and the power 
supply unit 4 so as to generate special pulse sequences which contain RF pulses and 
gradients. The phase and amplitude obtained from the demodulation unit 10 are applied to a 
processing unit 1 2. The processing unit 12 processes the presented signal values (also 
referred to as k-space) so as to form an image by transformation. This image can be 
visualized, for example by means of a monitor 13. 

Figure 4 is illustrative of a method for k-space data acquisition in accordance 
with the invention. The starting point for the k-space data acquisition is the central point 400 
of the k-space 402. From central point 400 the trajectory goes to the left in the -kx direction 
and then to the lower left comer of the k-space 402 in the -ky direction. 

Taking the lower left comer of the k-space 402 as a starting point a sampling 
404 is performed in the kx direction. After that the ky position is incremented by applying a 
corresponding gradient pulse. Further a compensation pulse such as a z-shimming pulse is 
applied. 

This way the trajectory continues in the k-space 406 where a sampling 408 is 
performed in the -kx direction. After the sampling 408 the ky poshion is further increased by 
applying a fiirther gradient pulse and another compensation pulse is applied such that the 
trajectory continues in the k-space 402. 

* 'In k-space 402 another sampling 404 is performed in the kx direction. After 

that again the ky position is incremented and a compensation pulse is applied. As a 
consequence the trajectory continues in the k-space 406 where another sampling 408 is 
performed in the -kx direction etc. By continuing this operation of alternately sampling the k- 
spaces 402 and 406 k-space data acquisition is performed for both k-spaces 402 and 406 
concurrently. 
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The resulting k-spaces 402 and 406 have both been imdersampled with 
alternating ky-lines being acquired. This resuUs in aliased images, which can be 'unfolded' 
using a multiple receiver coil set-up and known coil sensitivity prolBles, i.e. by applying 
parallel imaging, e.g. the SENSE or the SMASH method, or its hybrids and variants. 

Preferably the k-space data acquisition is performed only partially for the k- 
spaces 402 and 406 in order to reduce the read-out time. 

It is to be noted that the TE is about the same for both k-spaces 402 and 406 
when TB is defined as the time from the start of the data acquisition to flie time when the 
central point 400 of tl?e k-space 402 and central point 410 of the k-space ^06 are reached by 
the trajectory. This way the contrast resolution of the resulting image is increased. 

Rather than applying a gradient pulse after each horizontal scan to increase the 
ky position this can also be done only after every second scan. For example after sampling 
404 has been performed no gradient pulse is applied in order to increase the ky position. 
Rather only a compensation pulse is applied before sampling 408 is carried out on the same 
ky position. After sampling 408 the ky position is incremented by means of a corresponding 
gradient pulse. 

Further it is important to note that the gradient pulse to increase the ky position 
and the compensation pulse can be applied concurrently. 

Figure 5 shows a diagram of the gradient channel signals for the embodiment 
of Figmre 4 where the ky position is incremented after every horizontal scan. 

The signal Gx is illustrative of the horizontal gradient pulses which are applied 
for horizontal sampling along the k^ or -kx axis. The signal Gx has negative pulses el and 
positive pulses e2. The signal Gx is constituted by an alternating sequence of negative pulses 
el and positive pulses e2. For example a negative pulse el corresponds to one sampling into 
the kx direction (cf sampUng 404 of Figure 4) whereas a positive pulse e2 corresponds to a 
sampling into the opposite direction (c£ sampling 408 of Figure 4). 

The signal Gy represents the gradient pulses which are applied for 
incrementing the ky positions of the data samplings. A Gy pulse is appUed after each el and 
after each e2 pulse in the embodiment considered here. 

The signal Gz is illustrative of compensation pulses, such as z-shimming 
compensation pulses, which are applied during the data acquisition in order to "switch" 
between the two k-spaces (cf k-spaces 402 and 406 of Figure 4). A compensation pulse is 
applied after each signal el or e2. 
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It is important to note, that the k-space lines for the two images obtained are 
sampled in an interleaved manner, and that as a consequence^ all read-outs for one of the k- 
spaces have the same sign. This is advantageous in EPI reconstruction of the individual 
images (cf. steps 5 14 and 5 16) as no phase correction is needed. The final image is obtedned 
as a SSQ or MIP of the two images of the two k-spaces. For the image construction il is a 
particular advantage that the echo times are equal with only one echo spacing difference. 
Figure-6-is-illu5trative-of-a-corresponding flow-diagram; 



In step 500 data acquisition is performed into the kx direction of the first k- 

space. In step 502 the ky position for the sampling is incremented by applying a 
10 corresponding gradient pulse. In step 504 a z-shimming pulse is applied which brings the 

trajectory into the other k-space. 

In the other k-space a sampling is performed into the opposite -kx direction. In 

step 508 ky is again incremented and another z-shimming pulse is applied in step 510 which 

brings the trajectory back into the first k-space. 
15 In step 512 another sampling is performed into the kx direction. This procedure 

continues until a sufficient amount of data samples for the two k-spaces has been obtained for 

image generation. 

In step 514 a first image is generated based on the data samples acquired in the 
first k-space by sampling into the kx direction. Likewise in step 516 a second image is 
20 generated based on the data samples acquired in the other k-space by sampling into the 
opposite direction. 

In step 520 the two images of steps 514 and 516 are combined to provide the 
resulting image. 

In essence. Figures 5 and 6 are descriptive of a method which employs a 
25 parallel imaging reduction factor of 2, and 2 z-shim values: echoes for the first and second 
image are sampled alternately, and consecutively blips along Gy move to the next ky line 
while moving to the next k-space. As such, each individual k-space is undersampled by a 
factor of 2. The reconstruction of the corresponding aliased images is performed using the 
known coil sensitivities according to e.g. the SENSE or the SMASH methodology. The z- 
30 shim gradient is applied between the echoes for the first and second image, and rewound 
before acquiring the next ky position echo for the first image. 

Figure 7 is illustrative of an alternative method for a k-space data acquisition. 
The same reference numerals are used in Figure 7 as in Figure 4 in order to designate like 
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elements. In contrast to the embodiment of Figure 4 a ky gradient pulse is applied only after 
every second horizontal scan in the k-space 406. 

Figure 8 shows a diagram of the representative gradient chaimel signals for the 
preferred embodiment where the ky position is incremented after every second horizontal 
5 scan* Li the embodiment considered here a Gy pulse is applied after each pair of el and e2 
pulses. 

Figure 9 is illustrative of an embodiment where an interleaved k-space data 
acquisition scheme is used which involves a number of n = 4 k-spaces. Figure 10 shows the 
coirespondiriFg gradient channel signalsi The echoes for the first to the forth image are 
1 0 sampled altematingly, and consecutive blips along Gy move to the next ky line while moving 
to the next k-space. The z-shim gradient is applied between the echoes for the first, second, 
third and forth image, and rewound before acquiring the next ky position echo for the first 
image. This way a parallel imaging reduction factor r of r=n=4 is accomplished. 
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LIST OF REFERENCE NUMERALS 



magnetic resonance device 1 

magnet system 2 

5 gradient coil 3 

supply unit 4 

traxrsmittBrcbtl -5 — 

modulator 6 

• object 7 

1 0 transmission / reception circuit 9 

demodulation unit 1 0 

control unit 11 

processing unit 12 

monitor 13 

15 central point 100 

comer 102 

central point 400 

k-spacc 402 

sampling 404 

20 k-space 406 

sampling 408 

central point 410 
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CLAIMS: 



1. 



A method for k-space data acquisition for magnetic resonance imaging (MRI) 



characterized by using of at least first and second k-spaces for interleaved data acquisition, 
the at least first and second k-spaces covering substantially the same physical region. 



second coordinate axis, the method comprising: 

a) sampling into a first direction along the first coordinate axis, 

b) ^plying a first compensation pulse, 

c) sampling into a second direction along the first coordinate axis, the second direction 
being opposite to the first direction, 

d) applying a second compensation pulse, 

e) repetitively carrying out the steps a) to d), 

3. The method of claims 1 or 2, the first and second compensation pulses being 
z-shimming pulses. 

4. The method of claims 1, 2 or 3, further comprising incrementing a sampling 
position on the second coordinate axis after each step a) and / or after each step b). 

5. The method of anyone of the preceding claims 1 to 4, further comprising 
performing a partial k-space data acquisition by means of the sampling in steps a) and c). 

6. The method of anyone of the preceding claims 1 to 5 further comprising 
combining the sampled data of steps a) and c) to generate an image. 

7. The method of anyone of the preceding claims 1 to 6 further comprising 
generating a first image based on the data samples being acquired in the first direction, 
generating a second image based on the data samples acquired in the second direction and 
combining the first and second images into one image. 



2. 



The method of'claim 1 , the k-spaces having a first coordinate axis and a 
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8. The method of anyone of the proceeding claims 1 to 7, whereby a number of n 

k-spaces is used for the interleaved data acquisition, and further comprising the steps of: 

- applying a number of n-1 first compensation pulses of a first amplitude, 

- applying the second compensation pulse with a second amplitude, where the second 
amplitude is n-1 times the first amplitude. 



9. A magnetic resonance imaging (MRI) device comprising means for 
interleaved k-space data acqxusition in at least first and second k-spaces. 

10. The magnetic resonance imaging device of claim 9 comprising: 

- means (2, 3, 4, 5, 6, 9, 10) for k-space data acquisition, the k-spaccs having a first 
coordinate axis and a second coordinate axis, 

- a control unit (11) for generating of control signals for the means for k-space data 
acquisition, 

wherein the control signals cause the means for k-space data acquisition to 

a) sample into a first direction along the first coordinate axis, 

b) apply a fii*st compensation pulse, 

c) sample into a second direction along the first coordinate axis, the second direction 
being opposite of the first direction, 

d) apply a second compensation pulse, 

e) repetitively carry out the steps a) to d). 

11. A computer program product for k-space data acquisition for magnetic 
resonance imaging (MRI), the computer program product comprising program means to 
perform k-space data acquisition in at least first and second k-spaces in an interleaved way. 

1 2. The computer program products of claim 1 1 , the k-spaces having a first 
coordinate access, the program means being adapted to perform the stq^s of : 

a) sampling into a first direction along the first coordinate axis, 

b) applying a first compensation pulse, 

c) sampling into a second direction along the first coordinate axis, the second direction 
being opposite to the first direction, 

d) applying a second compensation pulse, 
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e) repetitively carrying out the steps a) to d). 



PHNL02094il||||^ 

14 01.10.2002 

ABSTRACT: 



The invention relates to a method for interleaved k-space data acquisition for 
magnetic resonance imaging (MRI), the k-spaces having a first coordinate axis and a second 
coordinate axis, the method coniprising: 

a) sampling into a first direction along the first coordinate axis, 

b) applying a first compensation pulse, 

c) sampling into a second direction along the first coordinate axis, the second direction 
being opposite to the first direction, applying a second compensation pulse, 

d) repetitively carrying out the steps a) to d). 



(Figure 4) 
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Sampling into kx direction 



Incrementing k% 



Apply z-shimming pulse 



Sampling into - kxdirection 



Incrementing kx 



Apply 2-shimming pulse 



Sampling into kx direction 



•500 



-502 



-504 



•506 



•508 



-510 



-512 



Generation of image on data samples acquired in kx direction 



Generation of image data on samples acquired in • kx direction 



Combination of images 



-514 



•516 



•520 



FIG. 6 
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k- space 2 



k • space 3 



k- space 4 
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